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Moving  Nonlocal  Crack:  Anti-Plane  Shear  Case 

N.  Ari+and  A.C.  Eringen 

Princeton  University 
Princeton,  NJ  08544 


ABSTRACT:  Stress  distribution  near  the  tip  of  a  constant  velocity 
crack  is  determined  by  means  of  the  nonlocal  theory  of 
elasticity,,  The  stress  at  the  crack  tip  is  finite  and 
it  increases  with  increasing  crack  velocity. 

I.  Introduction 

The  present  work  is  concerned  with  the  investigation  of  the 
stress  distribution  near  the  tip  of  a  uniformly  moving  crack  in  a  brittle 
elastic  solid.  To  this  end,  we  employ  the  recently  developed  theory  of 
nonlocal  elasticity  [1]  which  incorporates  important  features  of  atomic 
lattice  dynamics  relevant  to  the  study  of  microscopic  defects,  as  well  as 
macroscopic  phenomena  that  fall  within  the  domain  of  classical  elasticity 

in  the  long  wave-length  limit. 
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In  a  series  of  papers,  Erin gen  and  his  coworkers  [2-6]  treated 
various  problems  on  crack  tip  stresses,  dislocations  and  dispersive  wave 
propagation  by  means  of  the  linear  theory  of  nonlocal  elasticity.  We 
mention  here  briefly,  the  following  few  basic  results  in  support  of  the 
power  and  the  potential  of  nonlocal  theory:  In  contrast  to  the  predictions 
made  in  classical  elasticity: 

(i)  Plane  waves  are  dispersive  and  the  dispersion  curves 
approximate  those  of  the  atomic  lattice  dynamics  and 
the  phonon-dispersion  experiments  within  a  6%  error 
margin  in  the  entire  Bril louin  zone  [1-3]. 

(ii)  The  hoop  stress  is  finite  and  acquires  a  maximum  near 
the  crack  tip*  Consequently,  physically  meaningful 
maximum  stress  hypothesis  may  be  used  at  all  levels  as 
a  fracture  criterion  [4-8]. 

(iii)  Theoretical  strength  of  solids  estimated  by  means  of 

nonlocal  theory  is  in  excellent  agreement  with  the  pre¬ 
dictions  based  on  atomic  theory  [9,10]. 


In  dynamic  problems  of  crack  propagation,  classical  theory  also 

predicts  infinite  stresses  at  the  crack  tip.  Moreover,  the  critical 

1/2 

dynamic  stress  intensity  factors  (i.e.  K  =  r  (t  ,  * yZ )  as 
0  =  0)  turn  out  to  be  independent  of  the  crack  velocity,  identical  to 
their  static  values  [11,12]. 

The  problem  of  a  constant  velocity  crack  offers  a  convenient 
testing  ground  for  the  extension  of  the  treatment  of  nonlocal  crack  problems 
to  dynamic  cases.  Here,  we  study  the  motion  of  a  finite  length  crack  for 
the  anti-plane  shear  case  (Mode  III).  This  problem  represents  an  oversimpli¬ 
fied  picture  of  the  dynamic  fracture.  Nevertheless,  it  does  offer  us  an 
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insight  on  the  nature  of  the  nonlocal  stress  field  near  the  tip  of  a 
moving  crack.  It  is  found  that  the  nonlocal  shear  stress  at  the  crack  tip 
is  finite  and  it  increases  with  increasing  crack  velocity.  These  charac¬ 
teristics  of  nonlocal  stress  behavior  enables  us  to  extend  the  maximum 
stress  hypothesis  to  the  dynamic  cases. 

In  Section  2,  we  summarize  basic  equations  of  nonlocal  elasticity. 
In  Section  3,  an  analytical  (approximate)  solution  is  obtained  for  the 
dual  integral  equations  ensuing  from  the  mixed  boundary  value  problem.  The 
nonlocal  stress  field  is  given  in  Section  4,  and  compared  with  the  local 
and  the  static  stresses. 


II.  Nonlocal  Elasticity 

For  linear  homogeneous  and  isotropic  elastic  solids,  the  two 
dimensional  stress  constitutive  equations  are  given  by  [1] : 


t 


k£ 


a(|x'-x|,e)  ok£(x’)  ds(x') 


S 


(2.1) 


where  S  is  the  two-dimensional  plane  region  and  a]cjj(x,»'t)  is  the  local 
(classical)  stress  tensor  at  x'  which  depends  only  on  the  local  strain 
tensor  e,  £(x',t)  at  the  point  x1  ,  where 


a 

e 


U 

k£ 


X  err  6kl  *  2lJ  eM 

7  (uk,l  *  ui,k5 


(2.2) 


(2.3) 
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Here  is  the  displacement  vector,  X  and  y  are  the  Lame  constants,  and  an 
index  following  a  comma  represents  a  gradient,  e.g* 

9x^ 

The  nonlocal  kernel  a(|x’-x|>E)  characterizes  the  range  and  the  strength 
of  interatomic  interactions .  e  is  the  nonlocality  parameter.  It  corre¬ 
sponds  to  an  internal  characteristic  length  of  the  material,  hence  reflects 
the  microstructure  and  discreteness  of  the  body.  The  nonlocal  kernels  may  be 
determined  by  requiring  that  the  wave  dispersion  equations  derived  from  the 
nonlocal  equations  approximate  the  corresponding  dispersion  relations  which 
are  derived  by  means  of  lattice  dynamics  or  obtained  experimentally.  In 
addition,  they  must  fulfill  certain  consistency  conditions  so  that  in  the 
case  e  0  (the  continuum  limit) ,  the  nonlocal  equations  revert  to  the 
classical,  local  equations,  i.e. 

(i)  lim  a(| x'-x| ,e)  =  6(x'-x) 

e  -*■  0 

(ii)  f  a( | x'-x| )  dS(x')  =  1  (2.4) 

J  S  ~  ~ 

(iii)  a  vanishes  as  |x'-x|  ■+ 

In  [5],  we  introduced  a  two-dimensional  kernel  which  satisfies  the 
consistency  conditions  and  also  yields  dispersion  curves  closely  approxi¬ 
mating  those  of  a  two-dimensional  perfect  square  lattice.  (For  other  possi¬ 
bilities,  see  [3,9]  ) 

a(|x'-x|,e)  =  (2tt  e2)-1  K0(e_1[(x'-x)2  +(y*-y)2]%  (2‘5) 

- 

We  note  that  the  homogeneous  kernel  a(|x'-x|)  used  here  violates  the 
inhomogeneity  that  exists  in  an  atomically  thin  layer  near  the  crack  boundary. 
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(2.5)  has  the  additional  convenient  property  that 

(1  -  e2V2)a  =  6(|x»-x|)  (2.6) 

For  perfect  crystals,  the  dispersion  curve  can  be  matched  exactly  at 
the  ends  of  the  Brillouin  zone,  leading  to  the  determination  of  e  : 

e  =  (0.22  -  0.31)  a  (2.7) 

where  a  is  the  atomic  lattice  parameter. 

In  the  sequel,  we  will  use  (2.5)  as  the  nonlocal  kernel  and  e  as 
given  by  (2.7)  so  that  there  will  be  no  parameter  adjustments  for  the  spe¬ 
cific  problem  treated  here. 

Cauchy's  equations  for  the  linear  momentum  remains  valid: 

^^k  Pu£  (2*8) 

where  p  is  the  mass  density  and  a  superposed  dot  denotes  the  material  time 
derivative. 

Upon  applying  the  operator  (2.6)  to  (2.1)  and  using  (2.3)  and  (2.4), 
we  obtain  the  following  singularly  perturbed  partial  differential  equations: 

(1  -  e2V2)  p'q  =  (2.9) 

In  Section  3,  we  obtain  the  solution  of  Eqs.  (2.9)  and  (2.1)  for  the 
nonlocal  moving  shear  crack  problem. 
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III,  Anti-Plane  Moving  Crack 

A  line  crack  of  length  2£  is  assumed  to  propagate  with  a  constant 
velocity  V  ,  in  an  elastic  plate  in  the  X^  =  0-plane.  The  uniform  motion 
of  the  crack  is  maintained  by  an  anti-plane  shear  stress  Tq  (Fig.  1).  In 
the  moving  coordinate  system 

x  =  X  -  Vt  y  =  Y  (3.1) 


we  have 


U1  =  0  >  u2  =  0  ,  u3  =  w(x,y)  =  w(-x,y) 


a 

xz 


a 

yz 


(3.2) 


(3.2) 


and  the  field  equations  (2.9)  reduce  to 


Wc2)2 


r ,  2  ,9 

[1  -  e  ( — j  + 

9x 


(3.3) 


where  C2  =  (d/p)  denotes  the  phase  velocity  of  shear  waves. 

We  note  that  the  nonlocal  field  equations  revert  to  the  local  equations 
for  £  =  0  and  to  the  static  equations  for  V  =  0. 

By  utilizing  the  Fourier  transform  in  the  x-direction 


f(k,y) 


(2r)% 


00 


f(x,y)  exp(ikx)  dk 

—00 


(3.4) 


we  obtain  the  general  solution  of  (3.3)  for  y  >  0  as 
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w(k,y)  = 


2  k 

(-) 

V7T 


A(k)  exp(-k  y(ek)y  -  ikx)dk 


(3.5) 


where 


2  ,  2  „2  2  2  2W,  2  2, ,2, 2, 

Y  =  (c2~V  -  e  V  k  )/  (c2  -  e  V  k  ) 


(3.6) 


and  the  inversion  contour  is  chosen  in  such  a  way  that 


]R  y  >  0 


The  boundary  conditions  for  the  moving  crack  problem  in  the  new 
coordinate  system  are  equivalent  to  the  static  case.  Here  we  follow  [5] 
in  determining  the  self-consistent  nonlocal  boundary  conditions,,*  For  the 
uniform  shear  load  Tq  ,  they  are  given  by 


w(x,0)  =  U 

|x| 

>  i 

°yztx>0)  ’  -  T0> 

|x| 

<  i 

(3.7) 

w(x,y)  =  0 

as 

2  2 

(x  +  y  )  oo 

Insertion  of  (3.5)  into  (3„2)  and  (3.7)  yields  a  set  of  dual  integral 
equations 


where 


oo 

r 

K  D(k)  Y(£K)  cos  fc^  dK  =  1  ;  |£|  <1 

o 

no 

r 

D(k)  cos  dK  =  0  ;  |£|  >1 

J 

0 


(3.8) 


The  use  of  a  different  set  of  boundary  conditions  are  explored  in  other 
work  (e.g.  [10]). 
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K=  x/l  ,  K  =  k£  6  =  e/£ 

DM  =  (^(y0/T0£2)  A(k) 


(3.  9) 


The  standard  Ansatz 


D(k)  = 


t1/2  h(t)  J  (<t)  dt 


0 


leads  to  a  single  Fredholm  equation 


(3.10) 


h(t)  + 


K(t,n)  h(n)dn 


=  t1/2  : 


o  <  t  <  1 


(3.11) 


where 


K(t,n)  =  (tn) 


1/2 


K[Q(K)  -  1)]  J0(Kt)  J0(<n)  dK  (3.12) 


QGO  = 


[(c2 -K2)/(d2 -k2)]1/2  ;  K*(c,d) 

i  [(k2  -c2)/(d2 -k2)]1/2  ;  K£(c,d) 


d  =  c „/6v  , 


c2  =  d2  -  ]/S2 


(3.13) 


The  kernel  K(t,n)  can  be  expressed  in  a  more  compact  form  by  reducing  (3.13) 
to  an  integral  around  the  branch  cuts  by  contour  integration  (see  also  [9], 
p.  286),  i.e. 

d 


K(t,n)  =  (tn) 


1/2 


f 


kQ(k) 


✓  H^Oct)  JQ (id)} 

Ih^Vh)  J0(Kt)f 


dK 


;  t  >0 


;  t  <n 


(3.15) 


where  denote  Hankel  functions, 

v 
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In  order  to  derive  (3.14),  we  first  evaluate  the  auxiliary  integrals 
(t  >  nj 


i . 
3 


CtrD 


h 


z[G(z)-l]  H^(zt)  J0(zn)dz  ;  3  =  1,2 


C. 

3 


G(z)  = 


■  { 

V  0 


Q(z)  ;  zeC. 


Q(z)  ;  zee, 


(3.15) 


where  a  bar  denotes  the  complex  conjugate  and  z  =  p  +  iq  is  a  point  in  the 
complex  plane*  The  contours  and  C2  are  shown  in  Fig*  2»  There  are  no 

singularities  within  these  contours  (1^  =  0,1^  =  0)  and  contributions  from 
the  circular  arcs  vanish,  due  to  the  asymptotic  behavior  of  the  integrands 
on  each  Cj  0  As  the  next  step,  we  recall  that 


H^Cz)  +  Hq2)  (z) 


Vz) 

i  Y0(z) 


G(z)  =  G(-z) 

H^Ciqt)  J0(iqn)  =  -  hJ2)  (- iqt)  JQ(-iqn) 


(3.16) 


and  we  obtain 


h  +  J2 


K(t,n) 


(tq) 


1/2 


(tn) 


1/2 


d 


j  QOO  J0(Kt)  JQ(Kn)dK 
c  d 

[  Q(K)  i  YQ(Kt)  JQ(icn)dK 


(3.17) 

0 


C 

(3.14)  follows  easily  from  (3*17)°  For  rj  >  t  ,  we  merely  exchange  the  roles 
of  n  and  t  and  recover  (3.14). 
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The  new  form  of  K(t,n)  activates  us  for  the  trial  solution 

d 


h(t)  =  t 1/2 [g  + 


P(s)  JQ(st)ds] 


(3.18) 


Substituting  (3.18)  into  (3.11),  we  obtain 


AB  t1/2  *  t1/2 


B(s)  p(s)  JQ(st)  +  pt 


1/2 


q(k)  hJ13  (<)  JQ(xt)dX 


+  t 


1/2 


dK  X  Q(<)  JQ(xt) 


P (s )  G(x,s)dx  =  t1^2  (3.19) 


A  =  1  - 


TT 


/K2-c2\l/2  dx 

\77)  ~ 


“  §  "  (1  -  ^)1/2  =  s2  (3.20) 

C2 


B(s)  =  1  - 


2 

TT 


U  7 

/  2  2\% 

/x  -c  \  xdx 

[  ,2  2 ]2  2 

\d  -x  'x  -s 


o  ; 


c  <  s  <  d 


(3.21) 


G(x,s)  =  [x  hJ1;,(x)  jq(s)  -  s  jx(s)  H^i:)(x)]/(x2-s2)  (3.22) 


We  set  AB  =  1  and  p(d)  =0  as  an  integrability  condition  and  utilize  the 
asymptotic  expansions  of  the  Bessel  function  to  reduce  (3.19)  into  a  simpler 
form  for  an  approximate  solution: 
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.1/2 


3(^) 1/2  exp(i  (k 


37T/43)  Q  (<)  J0(Kt)  dK 


-  t 


1/2 


d 

f 


J 

c 


i  exp(i  (k-s)) 
^(KS)1/2  (<-S) 


p(s)  ds]  K  Q(K)  JQ(Kt)d< 


+  t 


1/2 


R(k3  Q(<) 


JQ(Kt)dK 


0 


(3.23) 


where  R(k)  is  given  by 


R(k) 


K 


d 

r 


[G(K.s)  ♦  l  3  p(s)  ds 

tt(ks)  ‘  (K-S) 

+  S  [H1  (k)  -  (^-)1/2  exp(i(K  -3tt/p)  )  ] 


(3.24) 


By  ignoring  the  last  term  in  (3.23),  we  obtain  an  airfoil  equation 
d 

—  p(s)  exp(-  is)  s  l/2  _  B(2/tt)^^2  exp(-  i 3tt/4) 

J  K  -  S  K 

c 


(3.25) 


The  solution  of  (3.25)  with  the  side  condition  p ( d)  =0  is  (Tricomi  [14], 
p.  178) 

p(s)  =  -  iB  (f)1/2  (^-)1/2  e=tp(i(s  -  31/4))  t|=|) 1/2  (3.26) 


Substitution  of  (3.26)  into  the  third  term  in  (3.24)  yields 
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d 

R(k)  Q(k)  J0(Kt)  =  r(t)  =  o(e1/2)  (3.27) 

c 

Hence,  p(s)  constitutes  a  very  good  approximate  solution  for  (3.11)  to  the 
order  of  (d^2)  . 

In  the  next  section,  we  utilize  the  approximation  provided  by  (3.18) 
and  calculate  the  nonlocal  crack  tip  stresses. 


IV,  Stress  Field 

The  nonlocal  stress  field  can  be  calculated  by  substituting  (3.18) 
into  (3.10)  and  (2.1).  Below,  for  comparison  purposes,  we  also  give  the 
crack  tip  stresses  for  the  local  and  the  static  cases.  The  stresses  along 
the  crack  line  near  the  crack  tip  (y  =  0  ,  £>1),  are  given  by: 

1 .  The  Static  Casej  V  =  0 

The  local  and  nonlocal  field  equations  coincide.  For  both  cases, 
the  solution  is  given  by 


D(k)  =  J1(k)/k 


(4.1) 


(i)  Local :  6=0 


LS, 

t  Tn 
yz  0 


J^(k)  cos  k£  dK  = 


I1(k)  e~K ^  dK 


a2  -  1)"1/2  [£  +  U2-l)1/2]_1  ;  5>1  (4.2) 


* 

The  approximate  solution  (5.26)  is  not  as  good  when  V  becomes  very  close 
to  c2  (i.e.  c  -K)) ,  since  then  the  asymptotic  expansions  of  the  Bessel 
functions  in  (3.23)  are  no  longer  valid. 
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The  second  integral  in  (4.2)  is  obtained  by  contour  integration  (i.e., 
J^k)  e1Ks  on  Cj  and  J^k)  e"1^^  on  C2  (see  Fig.  2)).  The  local 
stress  is  unbounded  at  the  crack  tip,  £  =  1. 

(ii)  Nonlocal :  6^0 


.NS  , 
yz  T0 


(1 +62k2)"%[(1  +62k2)4'  +  6k]'1  J2(<)  cosk£  d«  (4.3) 


.-1 


Ij(k)  e  Kc>  dK  + 


-1 


[1  -6k/(62k2-1)%]  I:(k)  e"K?  d< 

(4.4) 


X 


[2(5-Ur% 


CC-D/26 


[2/(£+l)]%  K  +  tt2-lfy~l 

X%  [K3/4(X)  *1/4  00  +  Kl/4 


-  1 

^  *3/4 


(X)]  +0(1) 

(4.5) 


From  (4.5)  and  from  the  small  argument  behavior  of  the  Bessel  functions, 
we  observe  that  the  crack  tip  stress  is  finite 

MC  1* 

tyz(1,°)/T0  =  +/2(rC3/4)/T(l/4))6£  =  0.475  (4.6) 

The  difference  between  (4.3)  and  (4.2)  lies  in  the  terms  which  include 
the  nonlocality  parameter  6  =  z/H  .  They  act  as  a  converging  factor  in 
(4.3)  and  the  nonlocal  stress  is  bounded.  Equation  (4.4)  is  obtained  by 
contour  integration  methods,  using  the  contours  and  C2  (Fig.  2). 

Here  we  see  how  the  nonlocality  acts  as  a  natural  cut-off  in  the  Fourier 
domain  and  thus  bringing  in  the  discreteness  effects  of  a  nonlocal  body. 
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2.  The Dunamia  Case:  V  /  0 


(i)  Local :  6=0 


The  solution  for  D(<)  is  different  from  the  static  case 


J,  OO 


D«)  =  (1  - 

C2 


_  x_r1/2  , 


J1(k)/(s2k) 


(4.7) 


Nevertheless,  the  local  stress  is  equal  to  its  static  value  and  it  is 
unbounded 


LD  . 

t  /x. 
yz  0 


s2kD(k)  cosk£  die  =  Jj(<)  cosk£  dK  (4.8) 


(ii)  Nonlocal :  6^0 


< 

ND  2  2  -ir  7  o  x.  -i 

t  /tq  =  3  (1+6  <  )  *  [(1+6  +  6k  Q(k)]"1  k  Q(k)  D(<)  cosk£  dK 


B(k)  Ij(k)  e'KS  d< 


*  (1-6k[(c^-V2  -  iW)/(c2  -  52vV)(iV  -1)]*)  B(K)  d< 


*  0(1)  ; 


where 


i  (i  -v2/c2)"1/4  +  (iiiK) 

*  2  ^C+IK'7 
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In  Fig.  3,  the  numerical  evaluation  of  (4.9)  are  plotted  along  with  the 
local  anti-plane  shear  stress.  In  contrast  to  the  local  stress,  the  non¬ 
local  stress  is  increasing  with  the  increasing  crack  velocity. 

In  this  work,  we  have  focused  on  how  the  nonlocality  affects  the 
stress  distribution.  In  an  upcoming  publication,  we  will  discuss  the  im¬ 
plications  of  the  velocity  dependence  of  the  nonlocal  crack  tip  stresses 
for  brittle  crack  propagation.  The  extension  of  the  maximum  stress 
hypothesis  to  the  dynamic  case  and  the  determination  of  a  terminal  velocity 
will  be  discussed  in  the  context  of  a  constant  velocity  crack  in  in-plane 
extension  mode,  which  provides  a  more  realistic  picture  of  dynamic  rupture 
phenomena. 


-i  o  o 
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“TV, 
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FIGURE  I 

MOVING  ANTI -PLANE  CRACK 


FIGURE  2 

CONTOURS  IN  THE  COMPLEX  PLANE 


FIGURE  3 

THE  NONLOCAL  STRESS 
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